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One of the great challenges in the ecology of infectious diseases is to understand what drives the emer- 
gence of new pathogens including the relationship between viruses and their hosts. In the case of the 
emergence of Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV), several studies have shown 
coronavirus diversity in bats as well as the existence of SARS-CoV infection in apparently healthy bats, 
suggesting that bats may be a crucial host in the genesis of this disease. To elucidate the biogeographic 
origin of SARS-CoV and investigate the role that bats played in its emergence, we amplified coronavirus 


wens F sequences from bat species captured throughout Thailand and assessed the phylogenetic relationships to 
Betacoronavirus : : ‘ fas 

Phylogeny each other and to other published coronavirus sequences. To this end, RdRp sequence of Coronavirinae 
Hipposideridae was targeted by RT-PCR in non-invasive samples from bats collected in Thailand. Two new coronaviruses 
SARS-CoV were detected in two bat species: one Betacoronavirus in Hipposideros larvatus and one Alphacoronavirus 
Thailand in Hipposideros armiger. Interestingly, these viruses from South-East Asia are related to those previously 
Emergence detected in Africa (Betacoronavirus-b) or in Europe (Alphacoronavirus & Betacoronavirus-b). These findings 


illuminate the origin and the evolutionary history of the SARS-CoV group found in bats by pushing for- 
ward the hypothesis of a Betacoronavirus spill-over from Hipposideridae to Rhinolophidae and then from 
Rhinolophidae to civets and Human. All reported Betacoronaviruses-b (SARS-CoV group) of Hipposider- 
idae and Rhinolophidae respectively cluster in two groups despite their broad geographic distribution 
and the sympatry of their hosts, which is in favor of an ancient and genetically independent evolution 
of Betacoronavirus-b clusters in these families. Moreover, despite its probable pathogenicity, we found 
that a Betacoronavirus-b can persistently infect a medium-sized hipposiderid bat colony. These findings 
illustrate the importance of the host phylogeny and the host/pathogen ecological interactions in the 
description and the understanding of pathogen emergence. The host’s phylogeny, biogeography and 
behaviour, combined with already described roles of pathogen plasticity and anthropic changes are likely 
to be co-factors of disease emergence. Elucidating the common ancestor of Hipposideridae and Rhinolo- 
phidae is key to understanding the evolutionary history of actual betacoronaviruses and therefore to get 
an insight of the deep origin of SARS-CoV. 

© 2011 Elsevier B.V. All rights reserved. 


1. Introduction could be the masked palm civet Paguma larvata (Viverridae) or an- 


other bridge-host, itself very probably contaminated by a SARS-like 


Pathogens mostly emerge in densely human populated areas as 
did the Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV), 
which had the potential to cause the first viral pandemic of the 
XXIst century (WHO, 2003; Jones et al., 2008). SARS-CoV most 
likely originated from bats through an intermediate species that 
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CoV from horseshoe bats (Song et al., 2005; Li et al., 2006; Cui et al., 
2007). Several elements indicate that this bridge-host could pre- 
sumably be a Carnivora as suggested by both detections in wild 
animals from markets or experimental infections (Guan et al., 
2003; Van den Brand et al., 2008). Several other mammalian or 
avian species are susceptible to coronaviruses infection (displaying 
enteric, hepatic, neurological, or respiratory tropism), but little is 
known about their viral ecology and host interactions. Coronavi- 
ruses are members of the order Nidovirales which contains three 
families of single stranded positive RNA viruses: Arteriviridae, 
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Coronaviridae and Roniviridae. Coronaviridae are divided in two 
subfamilies: Coronavirinae and Torovirinae. Because of serological 
and genotypic differences, coronaviruses were traditionally 


divided into three groups labeled group I, II and III (Gonzalez 
et al., 2003; Gorbalenya et al., 2004). These three groups have re- 
cently been elevated at the genus level and named respectively 
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Fig. 1. Phylogenetic reconstruction of main coronaviruses based on the analysis ofof 333 bases of the RNA-dependent RNA polymerase protein (partial nsp12) and 
geographical localization of Hipposideros larvatus Betacoronavirus Tata and Hipposideros larvatus Alphacoronavirus Ratcha. (A) Bayesian phylogeny of coronaviruses showing 
host taxonomy of Betacoronavirus-b group and reporting two new coronaviruses (posterior probabilities are shown for each node). (B) Repartition of hosts families (Hutson 
2003), sampling sites and HI BCoV T & Ha BCoV R spatial position. 
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Alphacoronavirus, Betacoronavirus and Gammacoronavirus by the 
coronavirus Study Group of the International Committee for 
Taxonomy of Viruses, thus invalidating the “Coronavirus” genus 
name (http://talk.ictvonline.org/media/p/1230.aspx). As other Nid- 
ovirales, Coronaviridae such as SARS-CoV possess a mutation prone 
RNA-dependent RNA-polymerase (RdRp) associated with high gen- 
ome mutability (Castro et al., 2005). Furthermore, they frequently 
recombine in vivo (Lai et al., 1985; Makino et al., 1986; Keck et al., 
1987; Jia et al., 1995; Herrewegh et al., 1998; Lee and Jackwood, 
2000; Stadler et al., 2003). This genetic variability confers to coro- 
naviruses a high potential of evolution that allows them to some- 
times overcome species barriers and host specificity (Baric et al., 
1995, 1997) (Thackray and Holmes, 2004). From a more general 
viewpoint, these characteristics are known to account for the fre- 
quent involvement of RNA viruses in disease emergence (Domingo 
and Holland, 1997; Chomel, 1998; Daszak et al., 2000; Cleaveland 
et al., 2001; Woolhouse and Gowtage-Sequeria, 2005). 

Until now, Asian SARS-related coronavirus (SARSr-CoV, from 
Betacoronavirus-b group —- Fig. 1A) RNA has been detected in only 
one bat family: the Rhinolophidae Gray, 1825 (Table 2). Therefore, 
this monogeneric bat family is considered as the host at the origin 
of the deadly SARS-CoV outbreak (Lau et al., 2005). However, other 
bat families are involved in the maintenance and the diversifica- 
tion of a wide range of coronavirus species belonging to both 
Alphacoronavirus and Betacoronavirus genera (i.e. Vespertillionidae, 
Phyllostomidae, Pteropodidae, Molossidae, and Hipposideridae 
families have been positively tested for coronavirus nucleic acid) 
(Chu et al., 2006; Tang et al., 2006; Dominguez et al., 2007; 
Carrington, 2008; Brandao et al., 2008; Tong et al., 2009; Pfefferle 
et al., 2009). Whether or not another bat family could be involved 
(in addition to Rhinolophidae) in the natural history of betacornav- 
iruses would have been an interesting question but since the dis- 
covery of SARSr-CoV in wildlife, most of the attention and 
sampling efforts have been focusing on Rhinolophidae (Table 2). 
Several species from this family were shown to carry different 
SARSr-CoV such as SARSr Rhinolophidae Bat Coronavirus strain Rp3 
in Rhinolophus pearsoni and SARSr Rhinolophidae Bat Coronavirus 
strain Rf1 in Rhinolophus ferrumequinum in the Chinese provinces 
of Guangxi and Hubei, respectively (Lau et al., 2005; Li et al., 
2005; Tang et al., 2006). The sampling pressure on both Rhinolo- 
phidae (1148 samples in 12 species - Table 2) and Vespertilionidae 
(900 samples in 17 species — Table 2) has induced an over repre- 
sentation of these families compared to others (Lau et al., 2005; 
Li et al., 2005; Tang et al., 2006). By contrast, Hipposideridae 
(206 samples in 6 species) were under-represented in previous 
studies (Poon et al., 2005; Woo et al., 2006). The higher northern 
distribution limit of Rhinolophidae (versus Hipposideridae) may 
partially explain the sampling bias in non-Southern Chinese prov- 
inces (Fig. 1B). This bias illustrates the need for testing other sym- 
patric species for betacoronaviruses to study their natural history 
and emergence process. 

The biogeography and the natural history of the Betacoronavi- 
rus-b species, among which are SARSr-CoVs, remain poorly de- 
scribed. However, Betacoronavirus have been detected in China, 
Africa, Europe and even Australia (Guan et al., 2003; Lau et al., 
2005; Li et al., 2005; Miiller et al., 2007; Tong et al., 2009; Quan 
et al., 2010; Drexler et al., 2010) (Fig. 1B). The biogeographical ori- 
gin of the ancestors of SARSr-CoVs (from Betacoronavirus-b sub- 
group) is therefore far more extended than the epicentre of the 
human epidemic (Southern China). Understanding the biogeogra- 
phy of this group would be essential for assessing the risk for 
emergence of another SARS-CoV like agent. South-China and its 
Pearl River delta combined various ecological and anthropic factors 
that have led the SARS-CoV to emerge but the distribution of the 
viral clade at the origin of SARS-CoV is much wider than this region 
and many locations still remain to be investigated. 


Emergence results from a combination of numerous factors. 
Environmental changes (ecological, climatic and agricultural), 
whether they are man-driven or not (Schrag and Wiener, 1995; 
Kuiken et al., 2003), globalization (trade or travel) and societal 
changes (land-use, demography), are among the most reported 
ones (Reperant, 2010), but others have been highlighted recently 
like the hosts’ phylogenetic relationships (Daszak, 2010; Streicker 
et al., 2010), the variation of biodiversity (Keesing et al., 2010) 
and species interactions in ecological communities whether be- 
tween hosts or between hosts and pathogens (Cleaveland et al., 
2001; Woolhouse and Gowtage-Sequeria, 2005; Reperant, 2010). 
Consequently, the biogeography and host-pathogen relationships 
of Betacoronavirus-b are pivotal to understanding their radiation. 
The diversity of the SARSr-CoV has been described by an increasing 
number of published studies and the number of discovered and se- 
quenced viral taxa has improved the accuracy of the phylogenetic 
reconstruction of this group. 

In this study we explored the biodiversity and the ecology of 
Coronavirinae, enriched the current phylogeny and challenged the 
biogeography of this clade. More particularly, we focused on the 
ecology of the Betacoronavirus-b sub-group (Fig. 1A) by testing 
hypotheses on the origination of these viruses not only in Rhinol- 
ophidae but also in their sister taxa Hipposideridae. Considering 
that Hipposideridae might host a diverse community of coronavi- 
ruses comparable to that of Rhinolophidae because of their phylo- 
genetic proximity and comparable diversity (Nunn et al., 2004; 
Simmons, 2007), it is very probable that the prevalence and diver- 
sity of coronaviruses hosted by this family has been underesti- 
mated because of biased sampling (Table 2 - Lau et al., 2005; Li 
et al., 2005; Tang et al., 2006). We hypothesize also that some 
bat taxa could play a refuge role for betacoronaviruses and there- 
fore constitute a source of infection for other species, consistently 
with the recent application of the theory of island biogeography 
(MacArthur and Wilson, 1967) to emerging pathogens (Reperant, 
2010). This study is particularly focused on, but not restricted to, 
Betacoronavirus in Hipposideridae in South-East Asia (Table 1). 


2. Materials and methods 


To test our hypotheses, we targeted, by RT-PCR, the Coronaviri- 
nae RdRp RNA sequence (more precisely a coding fragment of a 
subunit of the RdRp, the non-structural protein 12 - nsp12 -) in 
bat samples and compare it with the published sequences down- 
loaded from public online databases. The SRAS-CoV nsp12 and 
nsp13 exhibit the highest conservation profiles among proteins 
from all coronaviruses (Stadler et al., 2003). Moreover, as RdRps 
(nsp12) are involved in replication and transcription processes 
and therefore play a central role in the virion offspring genome 
synthesis (the vertical information transmission), its relatively 
conserved ribonucleic acid sequence is appropriate for both viral 
detection and phylogeny reconstruction. 


2.1. Sampling 


From 2006-2008, samples (feces, urine and oral swabs) were 
collected from bats at 26 sites in Thailand (Fig. 1B and Table 1). 
Sampling locations were selected considering land cover and geo- 
logical data as forests are known to enhance biodiversity (Harvey 
et al., 2006; Trisurat et al., 2010) and limestone hills harbour 
numerous caves. The rocky forest-covered hills surrounded by cul- 
tivated lowlands, considered as potential refuge areas, were partic- 
ularly sought. Sampling sessions covered an area of approximately 
800 km from North to South Thailand by 300 km from West to 
East. Dates and locations are described in Table 1 and Fig. 1B. Sites 
were sometimes sampled several times to confirm the results, to 


M. Ar Gouilh et al./ Infection, Genetics and Evolution 11 (2011) 1690-1702 


1693 


Table 1 
Detection of coronaviruses in saliva (Oral Swab), feces and urine in bats from Thailand. 
Date Family Genus Species O Swab Faeces Urine [Province Cave 
N + N +(%)| N + 
2008-01 Craseonycteridae Craseonycteris thonglongyai - - 3 ie} 5 0 |Khanchanaburi 27 
2008-01 Craseonycteris thonglongyai - - 7 ie} 6 0 |Khanchanaburi 28 
Craseonycteridae cumulation| - : 10 0 11 ) 
2006-03 Hipposideridae Hipposideros _ lekaguli 1 0 - - 2 0 |Nan 1 
2006-03 Hipposideros __larvatus 4 0 - - 3 0 |Nan 1 
2006-03 Hipposideros —_armiger 2 ie} - - 2 0 |Nan 1 
2006-04 Hipposideros —_armiger 6 0 7 ie} 3 0 |Phayao 2 
2006-04 Aselliscus stoliczkanus 1 0 1 ie} 1 0 |Phayao 2 
2006-04 Hipposideros __larvatus 2 ie} 2 0 1 0 |Phayao 2 
2006-04 Hipposideros __lylei 2 0 2 e) 2 0 |Phayao 2 
2006-04 Hipposideros __larvatus 15 0 18 0 5 0 |Chiang Rai 3 
2006-04 Aselliscus stoliczkanus 2 ie} 2 ie} - -  |Chiang Rai 3 
2006-06 Hipposideros —_armiger - - 17 ie} - -  |Ratchaburi 4 
2006-06 Hipposideros larvatus * “ 7 ie} 1 Q |Ratchaburi 4 
2007-04 Hipposideros —_armiger - - 5 ie} - -  |Uthai Thani 13 
2006-11 Hipposideros __larvatus 12 0 - - 3 0 |Ratchaburi 5 
2006-11 Hipposideros armiger 2 0 1 ie} : -  |Ratchaburi 6 
2006-11 Hipposideros _lekaguli 9 0 1 ie} 2 0 |Ratchaburi 6 
2007-05 Hipposideros __ larvatus - - 8 0 - |Uthai Thani 14 
2007-05 Hipposideros —pomona - - 1 ie} - -  |Uthai Thani 14 
2007-05 Hipposideros cineraceus . - 2 0 - - |Uthai Thani 14 
2007-05 Hipposideros __larvatus - - 3 ie} - - |Tak 15 
2007-05 Hipposideros —_armiger - - 5 0 - - |Tak 16 
2007-05 Hipposideros __lylei - - 2 0 - -  |Tak 16 
2007-09 Hipposideros —_armiger 10 oO | 10 gay - - |Ratchaburi 4 
2007-09 Hipposideros __larvatus 1 0 ll ie} - -  |Ratchaburi 4 
2007-03 Hipposideros —_armiger 10 ie) 10 0 : - |Chachoengsao 17 
2007-03 Hipposideros _larvatus 1 of] 11 Bey - - |Chachoengsao 17 
2007-11 Hipposideros — pomona - - 3 ie} - - |Nakon Ratchasima 18 
2007-12 Hipposideros _cineraceus - - - - 5 0 |Khanchanaburi 24/25 
2007-12 Hipposideros —pomona - - 2 0 1 Q |Khanchanaburi 24/25 
2007-12 Hipposideros _halophyllus - - 2 0 4 0 |Khanchanaburi 24/25 
2008-09 Hipposideros __larvatus - - | 40 Par - - |Chachoengsao 17 
Hipposideridae cumulation| 90 0 | 163 35 OO 
2006-03 Megadermatidae Megaderma lyra 2 0 - - 1 0 |Nan 1 
2006-04 Megaderma spasma 5 ie} 5 ie} 3 0 |jLampang 12 
Megadermatidae cumulation| 7 0 5 0 4 0 
2007-05 Pteropodidae Eonycteris spelaea - - 2 ie} - -  |Tak 15 
2007-11 Rousettus leschenaultii - - 35 0 11 0 |Lopburi 21 
2007-11 Eonycteris spelaea - - 1 ie} - - |Nakon Ratchasima 18 
Pteropodidae cumulation| - E 38 0 11 t) 
2006-03 Rhinolophidae Rhinolophus yunanensis 1 0 - - - -  |Nan 7 
2006-03 Rhinolophus lepidus 3 ie} - - 2 0 |Phitsanulok 8 
2006-04 Rhinoloplus malayanus 1 ie} 1 ie} 1 0 |Phayao 2 
2006-04 Rhinoloplus sp. 3 0 2 0 2 0 |Phayao 2 
2007-04 Rhinolophus coelophylius - - 1 0 - -  |Uthai Thani 13 
2006-11 Rhinolophus thomasi 1 ie} - - 1 0 |Ratchaburi 5 
2007-05 Rhinolophus coelophyllus - - 6 ie} - - |Uthai Thani 14 
2007-05 Rhinolophus malayanus - - 1 ie} - -  |Uthai Thani 14 
2007-05 Rhinolophus thomasi - - 1 ie} - -  |Uthai Thani 14 
2007-05 Rhinolophus coelophyllus - - 5 ie} - - |Tak 15 
2007-05 Rhinolophus stheno - - 1 ie} - - |Tak 15 
2007-11 Rhinolophus pearsoni - - 5 ie} - -  |Saraburi 19 
2007-11 Rhinolophus pearsoni - - 1 ie} 1 0 |Lopburi 20 
2007-12 Rhinolophus malayanus - * 1 (0) * -  |Khanchanaburi 23 
2007-12 Rhinolophus malayanus - - 1 0 - -  |Khanchanaburi 23 
2007-12 Rhinolophus malayanus - - 2 ie} 1 0 |Khanchanaburi 24/25 
Rhinolophidae cumulation| 9 Oo | 28 0 8 0 
2006-03 Emballonuridae Taphozous longimanus 25 0 - 7 12 0 |Lampang 9 
2006-03 Taphozous longimanus 12 0 - - 4 0 |Nan 10 
2006-03 Taphozous longimanus 9 0 - - 1 0 |Nan 7 
2006-04 Taphozous longimanus 22 ie} - - - - |Phayao 2 
2006-04 Taphozous longimanus 5 0 5 ie} 1 0 |Chiang Rai an 
2006-04 Taphozous longimanus - - 3 ie} 3 0 |Phayao 2 
2007-12 Taphozous longimanus 2 0 2 0 - - |Chachoengsao 22 
2007-03 Taphozous longimanus 4 0 4 ie} 4 0 |Chachoengsao 17 
Emballonuridae cumulation| 79 t) 14 0 25 ) 
2006-03 Vespertilionidae Myotis siligorensis 1 0 - - - - |Nan 7 
2006-04 Myotis siligorensis 1 0 - - - - |Phayao 2 
2006-11 Myotis siligorensis - - 4 0 - -  |Ratchaburi 5 
2006-04 Myotis siligorensis 1 0 1 ie} - 0 |Chiang Rai 3 
2008-01 Myotis siligorensis - - - - 1 QO |Khanchanaburi 26 
Vespertilionidae cumulation| 3 0 5 0 1 (0) 
2006-04 Miniopteridae Miniopterus schreibersii 2 ie} 2 ie} 2 0 |Chiang Rai 3 
Total number of samples 
All families cumulation 190 265 97 552 
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screen another population or species of bats or to assess a putative 
long-lasting infection of a given colony (Table 1). Samples were 
collected from 25 species of bats and were divided and preserved 
in two different media: (i) modified Minimum Essential Medium 
antibiotic supplemented (viral transportation medium) to preserve 
entire virus particles integrity and (ii) RNAlater (RNA preservation 
medium) to limit RNA degradation, both kept on ice before and 
immediately after collection (prior to liquid nitrogen conservation 
in the field). Each sample was immediately aliquoted into four 
tubes (including absolute backup) and kept at —80 °C before RNA 
extraction. Traceability was assured by assigning a unique number 
to each sample. Handling of animals was conducted in accordance 
with the American Society of Mammalogists Guidelines for the 
capture, Handling and Care of Mammals (Gannon and Sikes, 2007). 
Trapping sessions were performed using harp-traps, flip-nets and 
hand-nets. Bats were transferred to individual clean cotton bags 
immediately after being captured and kept in a dark, cool and quiet 
place for one to a few hours. Immediately after sample collection, 
bats were released where they were trapped. Particular attention 
was paid: (i) to avoid disturbance inside the cave and near the col- 
ony; (ii) to disinfect material and shoes between sampling sessions 
and study sites; and, (iii) to preserve the surroundings as first 
encountered. 


2.2. RNA extraction 


Qiamp viral RNA mini kit (Qiagen) was used for RNA extraction 
following the manufacturer’s instructions. RNA was eluted in 60 jl 
RNAse-free saline buffer, 7 1 were immediately used for the re- 
verse transcription reaction and the remaining 53 ul were ali- 
quoted and kept at —80°C for subsequent detection assays. The 
amount and quality of the extracted RNA was measured using 
Nanodrop and Agilent techniques prior to reverse transcription 
(Mueller et al., 2004; Gallagher and Desjardins, 2006). 


2.3. Reverse transcription (RT), amplification by polymerase chain 
reaction (PCR) and sequencing 


Complement DNA synthesis was performed with M-MLV 
Reverse Transcriptase (Bio-Active) and random hexamer primers. 
Random hexamer primers were used at a final ratio of 0.15 pg 
for 1 ug of RNA. Then, 5 jul of RT products were used in 25 ul final 
volume PCR reaction targeting a 438 nucleotide fragment of the 
RdRp coding region. PCR primers were designed from the RdRp 
coding region by selecting a highly conserved region (nsp12) 
among all Coronavirinae sequences available from online public 
databases in August 2007. Assuming a higher mutation rate for 
the third codon position and according to the published sequence 
variations, degenerate primers were designed: BatCoV pol 15197 
(forward primer: 5‘’-GGTTGGGAYTAYCCWAARTGTGA-3’) and Bat- 
CoV pol 15635 (reverse primer: 5’-CCATCRTCMGAHARAATCAT 
CATA-3’), respectively at position 15197 to 15219 and 15635 to 
15614 on Bat SARS-CoV Rf1 reference sequence (DQ412042). In or- 
der to increase both the sensitivity and the specificity of the detec- 
tion, a third primer were designed and combined with the previous 
ones in a semi-nested PCR: BatCoV pol nested 15419 (forward 
nested primer: 5’-GCNAATWSTGTNTTTAACAT-3’), at the position 
15419 to 15438 on the same reference sequence. The first PCR pro- 
duced a fragment of 438 nucleotides and the semi-nested PCR con- 
firmed the detection and enhanced significantly the detection limit 
of the semi-nested PCR by amplifying a 216 base pairs fragment 
from the first amplicon. Contaminations were avoided by observ- 
ing strict lab procedures, carrying out at each step (mix, PCRI, 
PCRII) in a UV cabinet and using bleached plastic-ware. Manipula- 
tions included extraction negative controls at the beginning and at 
the end of each sample series as well as several PCR negative con- 


trols per run. The sequencing of the 438 nucleotides PCR products 
was performed in both directions using the external primers. 


2.4. Cell culture and inoculation assay 


Fifty microliters of sample preserved in modified Minimum 
Essential Medium antibiotic supplemented (viral transportation 
medium) were inoculated into 12 wells plates where subconfluent 
cell monolayers were grown. Madin-Darby Canine Kidney (MDCK), 
Vero and Vero clone E6 cells were used in these assays. MDCK were 
chosen because Carnivora are highly suspected to have efficiently 
contributed to the transmission from bats to Humans (Paguma 
larvata, Nyctereutes procyonides, Melogale moschata) (Guan et al., 
2003) and several species are susceptible to SARS-CoV infection 
(Martina et al., 2003). The Vero cell line is known to efficiently rep- 
licate SARS-CoV (Qinfen et al., 2004). Cells were cultivated in Min- 
imum Essential Medium supplemented with 120 ,1g/ml of 
streptomycin, 120 units/ml of penicillin and 10% of Foetal Calf 
Serum (FCS). Once confluence was nearly reached, the supernatant 
was removed and samples were added on the cell monolayer and 
left at 20 °C ona rocking agitator for 60 min. Then, 2 ml of the same 
medium (2% FCS) were added to each well. Plates were then incu- 
bated at 37°C and daily checked on a phase contrast inverted 
microscope for CPE occurrence. 


2.5. Phylogenetic reconstruction 


A 304 or 333 nucleotide region was selected within the se- 
quences from both amplicons (Ha AcoV R for 304 or HI Bcov T, 
respectively). The corresponding region from 60 taxa available in 
databases were added to the matrix used in the phylogenetic 
reconstructions. 

1/ Firstly a dynamic homology approach implemented in Poy 
v.4.1.2 was used (Wheeler, 1996; Janies and Wheeler, 2002; Var6n 
et al., 2010). This method did not need any prior alignment because 
this step was performed during the cladogram diagnosis and there- 
fore established the putative homology statements during the tree 
reconstruction. Poy was used under parsimony as the optimality 
criterion (Wagner algorithm), complemented with the tree bisec- 
tion reconnection swapping algorithm (TBR), perturbation injec- 
tion, ratcheting and fuse iterations steps (Nixon, 1999). These 
combined steps provided a reasonably fast and accurate way to 
better explore the tree space by escaping from local optima. 
Although RdRp is a conserved gene and the matrix included a great 
number of taxa, resampling taxa tests were used to evaluate 
whether this matrix was subjected to the long branch attraction 
phenomenon which could constitute a limitation of some phyloge- 
netic reconstructions. This limitation could be particularly true in 
the case of high evolution rate and very limited number of taxa 
when combined with static homology approach of molecular char- 
acters. The over-representation of sub-groups by several similar 
sequences was avoided as it could artificially report changes on 
branches located between internal nodes, leading to the long 
branch attraction issue. 

2/ Secondly, a probabilistic approach (Bayesian) implemented 
in BEAST (Drummond and Rambaut, 2007) was applied to the same 
data set. Sequences were previously aligned using Sea-View (Gouy 
et al., 2010) under the Muscle algorithm (Edgar, 2004) and resulted 
in a frank border matrix showing a unique gap of three codons for 
four sequences (Bulbul, Trush, Munia and Asian leopard CoVs, po- 
sition 33 to 36). Then the adequacy of the matrix to 88 substitution 
models was evaluated using jModeltest (Posada, 2008). The three 
codons partition model of evolution was specified (under HYU sub- 
stitution model with a gamma site heterogeneity distribution), so 
that each of the three codon position could follow its own evolu- 
tion rate. Yule speciation model showed similar estimates as the 
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Table 2 

Asian bat species sampled in coronavirus studies. Some values should be considered cautiously because the related data were not always available and were in some cases 
extrapolated: (*) Details not available, supposedly two samples per bat (respiratory and fecal numbers extrapolated here). (*) Expressed as number of bats in original publication, 
details not available. 


Date Region Family Genus species PCR Detection Serology Viral ID & Taxonomy (partial) Reference 
All Respiratory Faecal 
N + % N + % N+ % N + % 


Pteropodidae 


2003 / 2004 HKSAR Cynopterus sphinx 30 15 0 15 0 ND Poon et al. 2005 
Mar 04 Nanning, Guangxi Rousettus leschenaulti 110 ND 110 0 84 | | 1,19 Li et al. 2005 
Maoming, 
Mar 04 Guangdong Rousettus leschenaulti 45 ND 45 0 42 0 Li et al. 2005 
Maoming, 
Mar 04 Guangdong Cynopterus sphinx 27 ND 27 ~(~O0 17 0 Li et al. 2005 
July 04 Nanning, Guangxi Rousettus leschenaulti 110 55 0 55 0 ND Li et al. 2005 
Dec 04 Nanning, Guangxi Rousettus leschenaulti 0 ND ND 58 | | 1,72 Li et al. 2005 
Apr 2004 / 
july 2005 HKSAR Cynopterus sphinx 4 2 0 2 0 ND Woo et al. 2006 
Apr 2004 / 
july 2005 HKSAR Rousettus lechenaulti 2 2 0 2 0 ND Woo et al. 2006 
Nov 04 / 
March 06 2 Chinese sites Cynopterus sphinx 12 6 0 6 0 ND Tang et al. 2006 
Nov 04/ 
March 06 1 Chinese site Rousettus leschenaulti 62 31.20 31.20 ND Tang et al, 2006 
Pteropodidae Total:| 396 111,00 293 0 201 2 1 
Hipposidpridae 
2003 / 2004 HKSAR Hipposideros armiger 8 4 0 4 0 Poon et al. 2005 
2003 / 2004 HKSAR Hipposideros pomona 6 3 40 3 0 Poon et al. 2005 
2003 / 2004 HKSAR Hipposideros armiger 24 12 0 12 0 ND 0 Lau et al. 2005 
Apr 2004 / 
july 2005 HKSAR Hipposideros armiger 26" 13 0 130 ND Woo et al. 2006 
Apr 2004 / 
july 2005 HKSAR Hipposideros pomona 36" 18 0 18 0 ND Woo et al. 2006 
Nov 04 / 
March 06 13 sites Hipposideros armiger 116 58 0 58 #0 ND Tang et al. 2006 
Nov 04 / 
March 06 1 site Hipposideros larvatus 6 3°40 3 0 ND Tang et al. 2006 
Nov 04 / 
March 06 2 sites Hipposideros pratti 18 9 #0 9 0 ND Tang et al. 2006 
Nov 04/ 
March 06 1 site Hipposideros pomona 2 1 0 | 0 ND Tang et al. 2006 
Nov 04/ 
March 06 2 sites Coelops frithi 12 6 0 6 0 ND Tang et al, 2006 
Nov 04 / 
March 06 1 site Aselliscus stoliczkanus 14 7 0 ta i*} ND Tang et al. 2006 
Hipposideridae Total:| 206 134 134 ND 


Miniopteridae 


2003 / 2004 HKSAR Miniopterus magnater 32 16 6,25 | 16 12,5 Poon et al. 2005 
2003 / 2004 HKSAR Miniopterus pusillus 38 19 26,32] 19 63,16 Poon et al. 2005 
Miniopterus schreibersii 
2003 / 2004 HKSAR Japanese 8 4 25 4 25 Poon et al. 2005 
2004 / 2005 HKSAR Miniopterus magnater 46 23 «0 23 «0 ND 0 Lau et al. 2005 
2004 / 2005 HKSAR Miniopterus pusillus 48 24.0 24 | 12,5 | ND 0 Co ee, et al. 2005 
Apr 2004 / 
july 2005 Site B HKSAR Miniopterus magnater 102* 51 0 51 1,96 | ND Woo et al. 2006 
Apr 2004 / 
july 2005 Site D HKSAR Miniopterus pusillus 50* 25 0 25 16 | ND Woo et al. 2006 
Apr 2004 / 
july 2005 Site E HKSAR ND Woo et al. 2006 
June 2005 Miniopterus magnater 25 0 | ND ND ND Chu et al. 2006 
Aug 2005 Miniopterus magnater 44 0 | ND ND Hcooromn assum et al. 2006 
Aug 2005 Miniopterus pusillus 1 0 ND ND Chu et al. 2006 
Dec 2005 Miniopterus magnater 25 0 ND ND Chu et al. 2006 
Dec 2005 Miniopterus pusillus n Oo | ND ND Chu et al. 2006 
Mar 2006 Miniopterus magnater 20 ND ND Chu et al. 2006 
Mar 2006 Miniopterus pusillus 10 Bo |ND ND Chu et al. 2006 
2005 15 sites (3) Miniopterus schreibersi 270 135 0 ND Tang et al. 2006 
Nov 04 Yichang, Hubei Miniopterus schreibersi 1 ND 10 Li et al. 2005 
Miniopteridae Total:| 579 20 345/297 7 2,36 | 298 40 13,42| 1 0 
Vespertilionidae 

‘Apr 20047 
july 2005 HKSAR Myotis chinensis 16 8 0 8 Oo ND Woo et al. 2006 
july 2005 H-HKSAR Myotis ricketti 46 23°20 23 4,35 | ND Woo et al. 2006 
Apr 2004 / 
july 2005 HKSAR Nyctalus noctula 14 7 0 3 0} ND Woo et al. 2006 
208 HZ |-H FUszeSstec Beacoenams 
july 2005 Sites B-C HKSAR Pipistrellus abramus 28 14 14 28,57 | ND Woo et al. 2006 
2003 / 2004 HKSAR Pipistrellus abramus 6 3°40 3 420 Poon et al. 2005 
2003 / 2004 HKSAR Myotis myotis 6 3420 3 0 Poon et al. 2005 
2003 / 2004 HKSAR Myotis ricketti 10 5 0 5 0 Poon et al. 2005 
2004 / 2005 HKSAR Myotis chinensis 3 0 Oo ND O Lau et al. 2005 
2004 / 2005 HKSAR Myotis ricketti 2 ie) 0 ND 0 Lau et al. 2005 
2004 / 2005 HKSAR Nyctalus noctula 2 0 ie) ND O Lau et al. 2005 
July 04 Tianjin Myotis ricketti 42 21 0 21 0 ND Li et al. 2005 
Nov 04 Yichang, Hubei Myotis altarium 1 ND 1 0 1°40 Li et al. 2005 
Nov 04 Yichang, Hubei Nyctalus plancyi 1 ND 1 0 1°40 Li et al. 2005 
2005 4 sites (1-Hainan) Pipistrellus pipistrellus 54 27 «(0 27 22,22) ND Tang et al. 2006 

8 sites (3-Ahnui, 
2005 Henan, Sichuan) Pipstrellus abramus 82 41.0 41 34,15) ND Tang et al. 2006 
2005 2 sites (1-Hainan) Scotophilus kuhlii 86 430 43 11,63} ND Tang et al. 2006 
Nov 04 / 
March 06 4 sites Myotis daubentonii 82 41.0 41 0 ND Tang et al. 2006 
Nov 04/ 
March 06 1 site Myotis mystacinus 2 1°40 1 0 ND Tang et al. 2006 

8 sites (4 Ahnui, 

Yunnan, 

Guangdong, 
2005 Guangxi, Jiangxi) Myotis ricketti 106 53 (0 53 24,53) ND Tang et al. 2006 


(continued on next page) 
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Nov 04/ 
March 06 2 sites Myotis chinensis 6 3 0 3 420 ND Tang et al, 2006 
Nov 04/ 
March 06 9 sites Myotis sp.a 160 80 0 80 0 ND Tang et al. 2006 
Nov 04/ 
March 06 2 sites Nyctalus aviator 12 6 Oo 6 0) ND Tang et al. 2006 
Nov 04/ 
March 06 3 sites Nyctalus noctula 34 17° (0 17° (0 ND Tang et al. 2006 
Nov 04 / 
March 06 1 site Scotomanes omatus 2 1 0 1 0 ND Tang et al. 2006 
Nov 04 / 
March 06 1 site Barbastella leucomelas 2 1 i?) | 0 ND Tang et al. 2006 
1 site (1- 
2005 Guangdong) Tylonycteris pachypus 28 14°00 14 [| 14,29) ND ee eer alam et al, 2006 
Nov 04 / 
March 06 1 site laio 16 8 0 8 i) ND Tang et al. 2006 
Nov 04/ 
March 06 1 site Murina leucogaster 10 5 0 5 0 ND Tang et al. 2006 
Apr 2004 / 
july 2005 Sites A, B HKSAR Tylonycteris pachypus 42 21 0 21 4 19,05 | ws 
Vespertillionidae Total:| 901 446 4 0,9 |448 49 10,94 
Rhinolophidae 
2003 / 2004 HKSAR Rhinolophus affinis 4 2 0 2 0 Poon et al. 2005 
2003 / 2004 HKSAR Rhinolophus pusillus 2 1 0 1 0 Poon et al. 2005 
2003 / 2004 HKSAR Rhinolophus rouxi 12 6 Oo 6 0 Poon et al. 2005 
2004 / 2005 HKSAR Rhinolophus affinis 4 2 0 2 0 ND 0 Lau et al. 2005 
2004 / 2005 HKSAR Rhinolophus sinicus 118 59 0 59 ie 44,07 | 18 le et al, 2005 
Nov 04 Yichang, Hubei Rhinolophus pusillus 15 ND 15 0 ND Li et al. 2005 
Nov 04 Yichang, Hubei Rhinolophus ferrumequinum 8 ND 8 125) 4 O Li et a/. 2005 
Nov 04 Yichang, Hubei Rhinolophus macrotis 11 3°40 8 12,5 | 7 71,43 Li et al. 2005 
Dec 04 Nanning, Guangxi Rhinolophus pearsoni 41 120 30 10 46 28,26 Li et al. 2005 
Dec 04 Nanning, Guangxi Rhinolophus pusillus 8 2 0 6 0 6 33,33 Li et a/. 2005 
Apr 2004 / 
july 2005 Rhinolophus affinis 14 10: 7 i} ND Woo et al. 2006 
Apr 2004 / 
july 2005 Sites D,J HKSAR Rhinolophus sinicus 236 118 0 118 El 19,49 | ND Woo eft al. 2006 
Apr 2004/ — Sites F, H, | 
july 2005 HKSAR ND Woo et al. 2006 
Nov 04/ 
March 06 27 sites Rhinolophus pusillus 101 101 0 101 0 ND Tang et al. 2006 
Nov 04 / 
March 06 2 sites Rhinolophus malayanus 30 15 0 15 0 ND Tang et al. 2006 
Nov 04 / 
March 06 9 sites Rhinolophus affinis 120 60 0 60 ND Tang et al. 2006 
Hee a. 
Henan, Sichuan, 
2005 Hubei, Shandong) Rhinolophus ferrumequinum 82 41.0 41 9,76 | ND Tang et al. 2006 
5 sites Rhinolophus thomasi 24 12 0 12 ND Tang et al. 2006 
16 sites (1- 
2006 Shandong) Rhinolophus sinicus 132 66 0 66 1,52 | ND Tang et al. 2006 
10 sites (1- 
2006 Shandong) Rhinolophus pearsoni 96 48 0 48 2,08 | ND Tang et al, 2006 
2004 11 sites (1-Hubei) Rhinolophus macrotis 76 38 «(0 38 2,63 | ND Tang et al. 2006 
Nov 04 / 
March 06 2 sites Rhinolophus rex 4 2 0 2 0 ND Tang et al, 2006 
Nov 04 / 
March 06 = 4 sites Rhinolophus luctus 8 4 0 4 0 ND Tang et al. 2006 
Nov 04/ 
March 06 1 site Rhinolophus osgoodi 2 10 1 0 ND Tang et al. 2006 
Rhinolophidae Total:| 1148 599 650 61 9,38 


constant population size coalescent model. Thus, the Yule specia- 
tion model was chosen as tree prior. The MCMC chain was set to 
10,000,000 states to get an adequate posterior ESS value superior 
to 300. 

The support of the nodes was evaluated via a bootstrap analysis 
(Felsenstein, 1985) for the direct optimization method and poster- 
ior probabilities we assessed for the Bayesian method (Fig. 1 and 
S1). Except for the sequences of Hipposideros larvatus Betacoronavi- 
rus Tata (HI BcoV T) and Hipposideros armiger Alphacoronavirus 
Ratcha (Ha AcoV R) which were newly identified in this study, all 
other sequences used in phylogenetic reconstructions were down- 
loaded from the National Center for Biotechnology Information 
database (http://www.ncbi.nlm.nih.gov/nucleotide/). Each Gen- 
bank sequence identification number is indicated directly on trees 
(Fig. 1 and S1). 


2.6. Rooting the trees 


In the case of coronavirus phylogeny, rooting with one member 
of the sister group (Bafinivirus or Torovirus, the two genus of 
Torovirinae sub-family) would have been a good choice to study 
the relation within the in-group Coronavirinae. Unfortunately, the 
region that we used here was too divergent and hardly comparable 
to Torovirinae sequence. However, as Torovirinae sequence can not 


be used to root the trees in this case, we used the results provided 
by a published study (Snijder et al., 2003), that used Torovirinae as 
root. This study, based on an amino acid matrix, inferred a Corona- 
viridae phylogeny, and verified the monophyly of Coronavirus and 
the basal branching of Gammacoronavirus clade. According to these 
conclusions, we chose an avian Coronavirus taxa to root our trees 
and therefore study the Alphacoronavirus and Betacoronavirus gen- 
era phylogeny. An additional taxon, from the same clade of the taxa 
used as root (Gammacoronavirus - Avian Coronavirus Beaudette), 
was included in the analysis to allow the polarization of the char- 
acters (ancestral versus apomorphic) along the tree. 


3. Theory 


3.1. Identifying the phylogenetic signal of filiation in Betacoronavirus 
metagenome 


A genome comprises ribonucleic acid fragments that can inde- 
pendently be submitted to various selection pressures of different 
intensity. The history of these fragments can be deeply shaped by 
the function they encode, eventually resulting in incongruent phy- 
logenetic signals. The only way to fully recover and clearly exploit 
the phylogenetic signal of a given fragment (by comparison to the 
putative homolog fragment of other taxa) is to analyze it 
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Fig. 2. Colonies of Hipposideros at roost, showing the scattered distribution of individuals, with a zoom on the face of one specimen (bar is equal to 1 cm). (A) Hipposideros 


larvatus (B) Hipposideros armiger. 


independently. Therefore, the combination of several fragments in 
a same analysis should be carefully thought and documented as 
mixing the information they provide can produce a blurred signal 
and therefore can lead to an inaccurate phylogeny. Here we fo- 
cused on a fragment of the polymerase sequence only for phyloge- 
netic reconstruction of coronaviruses because: (i) the key function 
of this protein involved in vertical transmission and consequently 
in viruses filiation is likely to regulate variation and provide a phy- 
logenetic signal that reflects the filiation and, (ii) its relatively short 
length, genome position and associated transcription characteris- 
tics make it unlikely to be a putative recombination site. Moreover, 
the analysis of this region by the RDP3 package (Martin et al., 2010) 
did not detect any putative recombination sites. 

During the inference of a species tree, the identification of the 
appropriate vertical transmission of information, is one of the 


classical pitfalls of phylogenetic reconstruction. Therefore, careful 
attention must be paid in the choice of a genetic fragment that re- 
flects the species phylogeny and not only the gene history. A given 
genetic fragment, coding for a positively selected function which 
appeared in one species, may benefit another one by horizontal 
transfer (recombination), increasing the fitness of the recombinant 
virus and its progeny. This phenomenon have happened also be- 
tween viruses that are actually recognized as members of 
separated families as described for the two independent recombi- 
nation events that occurred between a Torovirus (Coronaviridae, 
Torovirinae), a Coronavirus and an Orthomyxoviridae (Influenzavirus 
C) (Snijder et al., 1991). These particular events did not involve the 
RdRp sequence. All these reasons lead us to rely on the analysis of 
this RdRp fragment for the phylogenetic reconstruction of 
Coronavirus. 
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3.2. Biodiversity, host phylogeny and micro hot-spots 


According to the distribution of both the known mammalian 
host families, the Rhinolophidae Gray, 1825, and its sister family 
the Hipposideridae Lydekker, 1891, Betacoronavirus-b might be 
widespread over Asia and most of the old world. Moreover, the bio- 
diversity of these two mammalian families is surely higher than 
what has been currently described (SJ. Puechmaille, unpublished 
data), which in turn induces a potential greater coronavirus diver- 
sity than is currently described. In discordance to what has been 
reported to date and based on the phylogenetic proximity of these 
sister families, Hipposideridae might harbour an equal or richer 
diversity of coronaviruses (at least of the Betacoronavirus genus) 
compared to the Rhinolophidae. Furthermore, the fragmentation 
of the forested habitat, produced by human activities such as land 
use practices and agriculture, leads to the isolation of preserved 
biodiversity spots (surrounded by cultivated fields) that have the 
potential to maintain both host and viral populations, preserving 
by this way the footprints of the natural history of coronaviruses 
in wildlife. 


4. Results 


Twenty-five species of bats belonging to the Rhinolophidae 
(N = 8), Hipposideridae (N = 8), Emballonuridae (N = 1), Molossidae 
(N = 1), Megadermatidae (N = 2), Pteropodidae (N = 2), Miniopteri- 
dae (N = 1) (Table 1), Craseonycteridae (N = 1) and Vespertilionidae 
(N= 1) families were sampled. Species were identified using mor- 
phological and acoustic characters. 

A total of seven samples among 552 (265 feces, 97 urine collec- 
tions and 190 oral swabs), were positive using the RT-PCR on the 
RdRp region described in the material and methods section. This 
amount rose to 28 when the results from the second sampling at 
the study site 17 were included (Table 1 and Fig. 1B). All of the po- 
sitive samples came from two species: Hipposideros larvatus 
(Fig. 2A) and Hipposideros armiger (Fig. 2B). Two new viruses were 
detected in this study but the isolation trials on cell cultures re- 
mained unsuccessful: no cytopathic effect was observed after three 
passages on cell cultures and no viral material was amplified from 
the third passage supernatants with the semi-nested PCR de- 
scribed in material and methods. 


4.1. Two coronaviruses from two Hipposideridae: Hipposideros 
armiger Alphacoronavirus Ratcha (Ha AcoV R) and Hipposideros 
larvatus Betacoronavirus Tata (HI BcoV T) 


The Ha ACoV R was detected in bat feces from a unique colony 
(Hipposideros armiger Hodgson, 1835). Among 10 bats tested, four 
were positive (36% - study site 4 —- Table 1). The BLAST 78% identity 
value indicated that the sequence was related to a bat Alphacorona- 
virus reported to infect Rhinolophus ferrumequinum and Rhinolo- 
phus pearsoni in the Chinese Shandong region in 2005 (Table 2) 
and (Tang et al., 2006). For this sequence, 30 bases from the 5’ side 
of the PCR products were ambiguous and therefore removed from 
the analysis because the sequencing with one degenerated primer 
gave low quality results. 

Concerning HI BcoV T, the sequences were free of ambiguous 
sites and identical among the molecular isolates amplified from 
feces collected from 22 specimens of bats from the same colony 
(Table 1). The 391 nucleotides long sequences obtained showed 
that this new virus was related to the SARS-CoV clade and shared 
75% nucleotide identity (online BLAST analysis) with the Bat SARSr 
Coronavirus Rf1 from Rhinolophus ferrumequinum. HI] BcoV T was 
detected in fecal samples from one colony of Hipposideros larvatus 
(Hipposideridae Lydekker, 1891), in one sampling site among 26. 


In this colony (estimated to 250 individuals), 3/11 (27%) bats tested 
were positive, then, eighteen months later, 21/40 (52%). This is the 
first evidence of SARSr-CoV long-lasting infection in a bat colony 
(Table 1). 


4.2. Other species and study sites 


Besides the detection of Ha ACoV-R in one site, no other Alpha- 
coronavirus was detected from the other study sites. Hipposideros 
armiger individuals positive for Ha ACoV-R shared their roost with 
Hipposideros larvatus individuals, which were negative for the 
detection of Ha ACoV-R. Similarly, besides Hl BcoV-T detected in 
one site, no other Betacoronavirus-b was detected from the other 
study sites even if they were populated by the same species as 
the positive site and in two instances (cave 3 and 14 - Table 1), 
with comparable sample sizes. The other species sampled, Taphoz- 
ous longimanus and Hipposideros armiger, that were sharing the 
same roost with Hipposideros larvatus (cave 17 - Table 1), were 
all negative for the detection of HI BcoV T. 


4.3. Phylogenetic reconstruction 


The same global topology was obtained from both direct opti- 
mization (Poy) and probabilistic methods (Beast) as apparent dif- 
ferences mostly concerned nodes that are not supported by 
respective bootstrap or posterior probability values (inferior to 
95%) (Fig. 1A and S1). Indeed, direct optimization supported Beta- 
coronavirus paraphyly but basal nodes exhibited low bootstrap 
support values. Accordingly, the separation between Alphacorona- 
virus and Betacoronavirus is not significantly supported in the 
Bayesian phylogeny (posterior probablility value: 0.69 < 0.95). Sev- 
eral other nodes remained weakly supported by the data in each 
method. Moreover, the HKU4/HKU5/BtCoV 133 group significantly 
segregated in basal position of the Betacoronavirus in both recon- 
structions. Alphacoronavirus exhibited mainly nodes with low sup- 
port in both methods of reconstruction despite comparable 
topologies. The Alphacoronavirus detected in this study, Ha AcoV 
R, rooted two alphacoronaviruses reported to infect bats in 
Bulgaria (Drexler et al., 2010). Concerning Betacoronavirus-b topol- 
ogy, most of the nodes were well supported. Our RdRp gene based 
phylogenetic reconstructions placed HI BCoV-T as a sister taxon of 
an african Hipposideros Betacoronavirus-b (Fig. 1A and S1). These 
taxa formed together the sister group of all other betacoronavirus- 
es-b, which comprised the monophyletic group of Rhinolophidae 
betacoronaviruses-b. These were rooted by one Hipposideridae 
Betacoronavirus-b: the Australian Bat Coronavirus R. aur, detected 
in Rhinonicteris aurantia. Even paraphyletic, Hipposideridae 
betacoronaviruses mainly grouped together and rooted the Rhinol- 
ophidae betacoronaviruses among which is the SARS-CoV (Fig. 1A 
and S1). 


5. Discussion 


The species demarcation criteria proposed to the ICTV by the 
coronavirus study group (http://talk.ictvonline.org/media/p/ 
1230.aspx) recommend a minimal 90% amino acid sequence iden- 
tity threshold over seven domains of the replicase (nsps 3, 5, 12, 
13, 14, 15, 16). Only one domain is available in our data (nsp12) 
but it is located in one of the three most conserved and contiguous 
regions of the replicase (nsp11, 12 and 13) (Stadler et al., 2003). 
The relatively low identity values obtained by the BLAST analyses 
(78% for the Ha AcoV-R, and 75% for the HI BcoV-T) in the most 
conserved region of the coronavirus genome suggests that HI BcoV 
T and Ha AcoV R are new species, even though there was only one 
fragment available, instead of the ad hoc set of seven regions pro- 
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posed by the ICTV working group. More investigations are 
necessary to fully characterize these coronaviruses and to get the 
identity values of the seven replicase domains to rule on these 
two taxa description in complete accordance with the ICTV recom- 
mendations. For clarity and non-ambiguous taxonomy, these two 
putative new species are proposed to be named Hipposideros larv- 
atus Betacoronavirus Tata (HI BcoV-T) and Hipposideros armiger 
Alphacoronavirus Ratcha (Ha AcoV R). These Betacoronavirus and 
Alphacoronavirus were detected from the same bat genus Hipposid- 
eros (super-family: Rhinolophoidea sensu Teeling et al. (2005), 
family: Hipposideridae Lydekker, 1891). Considering the diversity 
of both Chiroptera and Coronavirinae, we explicitly implied the 
host species and the virus species within the name, thus avoiding 
further possible confusion once other coronaviruses will be discov- 
ered and described within these genera (Betacoronavirus and 
Hipposideros). 


5.1. Sampling pressure and host species involved 


Sampling bias could explain the uneven distribution of corona- 
virus detection among bat families, more specifically concerning 
SARSr-CoV in Rhinolophidae versus Hipposideridae. In contrast 
with other studies, the Hipposideridae family in our study was 
the most sampled family with five times more samples than those 
collected from Rhinolophidae. Only two sampling sites among 28 
(site 17 for HI BcoV-T and site 4 for Ha AcoV-R) explored with var- 
ious sampling intensity, gave positive results. For some site, the 
low level of sampling might be insufficient to detect any Coronavir- 
inae (Table 1). We cannot exclude that betacoronaviruses are 
hosted by others hipposiderid bat colonies and could circulate 
around Thailand or South-East Asia. Moreover, it is likely that other 
yet unknown betacoronaviruses circulate in Hipposideridae and 
even in other bat families around South-East Asia. An appropriate 
sampling would be able to test these hypotheses. 


5.2. Long-lasting infection in colony, roost-sharing behaviour and host 
phylogenetic affinities are suspected to be species-jumping promoting 
factors 


Chu and collaborators have repeatedly detected alphacoronavi- 
ruses in Miniopterus (Miniopteridae) colonies from Hong Kong. 
Positive detections in bat colonies were obtained over a period 
ranging from 3 months in a Miniopterus pusillus colony to 6 months 
in a Miniopterus magnater colony (Chu et al., 2006). These results, 
suggested that alphacoronaviruses were able to persist through 
time in bat colonies. Another study on Alphacoronavirus (and Beta- 
coronavirus — SARSr-Rh-BatCoV —) infected Rhinolophus sinicus bats 
modulated this conclusion. No significant persistence was ob- 
served in Rhinolophus sinicus bat colonies and the viral clearance 
was estimated between 2 weeks and 4 months (Lau et al., 2010). 
Moreover, the authors reported a significant loss of body weight 
in bats infected by Betacoronavirus (SARSr-CoVs) that was not ob- 
served in Alphacoronavirus infected bats. This loss of weight could 
be the sign of a loss of fitness associated with Betacoronavirus 
infection in Rhinolophus sinicus, without presuming that it would 
be a promoting factor or a consequence of Betacoronavirus infec- 
tion. Even though: (i) Betacoronavirus persistence was not observed 
in Rhinolophus sinicus colonies and, (ii) the loss of body mass corre- 
lated with Betacoronavirus infection might reduce the fitness of 
Rhinolophus, we report here that a Betacoronavirus was persisting 
for 18 months in a Hipposideridae bat colony. This suggests that 
a Hipposideridae colony might be more tolerant to Betacoronavirus 
than Rhinolophidae over a long period of time or that Betacorona- 
virus of Rhinolophidae might have acquired virulence factor limit- 
ing long-lasting infection in this host. Further investigations are 
needed to confirm and describe this suspected difference. This 


Betacoronavirus detected so far in Hipposideros larvatus in Thailand 
is restricted to a single cave. Thus, Betacoronavirus long-lasting 
infection in a bat colony might be a rare event. Although rare, 
long-lasting infection can promote cross-species transmission by 
the emergence of new properties in viral populations (Baric et al., 
1999). Despite the coincidence of several promoting factors in 
study site 17 (the long-lasting infection - Table 1 -, the sharing 
of the same small-sized cave by three bat species and the phyloge- 
netic proximity of the two Hipposideros species) none of the other 
species were detected positive for Betacoronavirus (Hipposideros 
armiger and Taphozous longimanus were negative — Table 1). Inter- 
estingly, we noted that the Taphozous longimanus (negative) colony 
was roosting separately from the Hipposideros larvatus (positive) 
cluster whereas the Hipposideros armiger (negative) colony occu- 
pied also another specific area of the cave. Similarly, in site 4 (Ta- 
ble 1) the Hipposideros armiger colony, infected by ACoV-Ha, was 
roosting separately from the other (negative) species. Again, this 
typical spatial segregation within a roosting site (e.g. cave) limits 
direct contacts between species and is suspected to account for a 
behavioral barrier to cross-species transmission in the wild, as pre- 
viously described for Lentivirus and Felidae (Troyer et al., 2008; 
VandeWoude et al., 2010). Nevertheless, we mentioned these 
observations here but these mechanisms still remain to be pre- 
cisely described in the complex relationship between Rhinolophoi- 
dea and betacoronaviruses. 

By contrast, previous studies of Asian coronaviruses that fo- 
cused on the species level, suggested that there was no strict asso- 
ciation between bat species and Alphacoronavirus species (Chu 
et al., 2006; Tang et al., 2006), and other studies on Betacoronavirus 
in Africa reached similar conclusions (Quan et al., 2010). Species- 
barrier crossing has been reported occasionally for Betacoronavirus 
in Asian, European and African Rhinolophidae or Hipposideridae 
but more frequently in Miniopteridae (Chu et al., 2006; Tang 
et al., 2006; Tong et al., 2009; Drexler et al., 2010). Moreover, Mini- 
opterus species exhibit an aggregative behaviour, and mixed spe- 
cies colonies are commonly observed as described in literature 
(Benda et al., 2003). Mixed species roosting, inducing direct contact 
between species and therefore promoting transmission, in synergy 
with viral persistence in a colony and the phylogenetic proximity 
of the hosts, is likely to promote cross species transmission of 
alphacoronaviruses. As a result, a non-specific host-virus associa- 
tion for Alphacoronavirus in Miniopterus follows the same model 
as what has been described for these bats and their ectoparasites 
(Bruyndonckx et al., 2009). The apparent higher frequency of 
cross-species transmission (or the absence of strict virus-host asso- 
ciation) in Miniopteridae compared to Hipposideridae and Rhinol- 
phidae might partially be explained by their different roosting 
behaviour. 


5.3. Phylogeny and biogeography: Relation between betacoronaviruses 
filiation and host distribution 


The two methods of phylogenetic reconstructions employed 
here supported comparable relationships within Coronavirinae 
(Fig. 1 and S1). The uncertainty associated with some relationships 
was often similarly reported by each method. These phylogenetic 
issues might be reduced by increasing the size of the analysed frag- 
ment, that would increase the number of informative characters 
and therefore provide more information to support nodes. Adding 
new taxa would also improve the phylogeny by providing probable 
currently lacking information, within or around groups that are not 
yet well supported: (i) Betacoronavirus sub-groups relationships 
might be improved by looking for new taxa related to the group 
HKU9/Bat CoV Philippine Dilim (weak posterior probability: 0.59 
and bootstrap value: 0.24 are reported for this node); (ii) Other se- 
quences from new Betacoronavirus (or taxa from a new group) 
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might help to clarify the relationship between Betacoronavirus, 
Alphacoronavirus and putative new groups as discordance is ob- 
served between phylogenies and none of the method provide sig- 
nificant support to the nodes that link Betacoronavirus and 
Alphacoronavirus (posterior probability: 0.69 and bootstrap values: 
0.21; 0.06; 0.24). 

Interestingly, our study revealed that the Betacoronavirus de- 
tected in Hipposideridae in South-East Asia was related to those 
that currently circulate in Africa, according to the analysis of a 
RdRp fragment. Our analyses strengthen the diversity of Betacoro- 
navirus-b clade and extend their broad distribution. Indeed, SARS- 
CoV emerged from this Betacoronavirus-b clade (Fig. 1A) which is 
present in Africa, Europe, Asia and even Australia. 

Quan and collaborators pushed forward the possible occurrence 
of migration events from Africa to Asia to explain the basal phylo- 
genetic position of an African Betacoronavirus (compared to Asian 
rhinolophidae SARSr-CoV) together with the apparent absence of 
SARSr-CoV in Asian Hipposideridae (Quan, 2010). With the detec- 
tion of Betacoronavirus in hipposiderids in Asia, this hypothesis, 
while still possible, becomes less necessary to explain phylogenetic 
and geographic patterns of Hipposideridae and Rhinolophidae 
betacoronaviruses. The clustering of Betacoronavirus-b that corre- 
lates with the host phylogeny rather than with virus geolocaliza- 
tion suggests that the phylogenetic relationships between 
betacoronaviruses-b is mainly driven by the host phylogeny 
(Fig. 1 and S1). Accordingly, the phylogenetic relationships of the 
HI BCoV detected here (from South-East Asia) with the SARSr- 
CoV (from Africa) and other betacoronaviruses support the follow- 
ing hypothesis. Betacoronaviruses-b ancestors, meaning 
SARSr-CoVs ancestors, could have been historically hosted by the 
common ancestor of the Rhinolophidae and Hipposideridae and 
could have later evolved independently in the lineages leading to- 
wards Rhinolophidae and Hipposideridae betacoronaviruses. This 
point of view is supported by the topology of Betacoronavirus clade 
that exhibits the monophyly of Rhinolophidae betacoronaviruses 
despite their broad geographic distribution (from Africa and cen- 
tral Europe to South-East Asia and Northern Australia - Fig. 1) 
and shows the close relationship of Hipposideridae betacoronavi- 
ruses. Nevertheless, it is unlikely that putative host switching be- 
tween Rhinolophidae and Hipposideridae betacoronaviruses 
(ancestors), eventually associated with ancient migration events, 
would still be detectable. More investigations are necessary to ad- 
dress these points and test corresponding hypotheses. As well, new 
investigations would be necessary to confirm this suggested pre- 
dominant influence of the host phylogeny on Betacoronavirus-b 
clustering and whether other lineage would exist in Hipposideri- 
dae, as suggested by the report of a Betacoronavirus from Australia 
(C. Smith pers. com. - Fig. 1). 

The position of the Hipposideridae Betacoronavirus from 
Australia (basal to the Rhinolophidae betacoronaviruses), suggests 
the existence of an intermediary subgroup (Fig. 1 and S1) between 
Hipposideridae and Rhinolophidae betacoronaviruses that sur- 
vived within a rare endemic Australian bat species. The phyloge- 
netic status of this host, Rhinonicteris aurantia, is not well defined 
but this species is described as intermediary between Hipposider- 
idae and Rhinolophidae, while considered an Hipposideridae 
(Hand and Kirsch, 1998; Simmons, 2007; Armstrong, 2006; Arm- 
strong com. pers.). Albeit this putative intermediary betacoronavi- 
rus group count a unique member to date and its associated 
posterior probability is weak, its phylogenetic position correlates 
with that of its host. Therefore, we are looking forward new inves- 
tigations that would address these elements and screen betacoro- 
naviruses in other Hipposideridae and Rhinolophidae. 

The ecological niche of these bats and their behaviour provide 
them with frequent opportunities to spread viral variants among 
other nocturnal mammals such as other Chiroptera or Carnivora 


Viverridae (civet family) (Song et al., 2005; Guan et al., 2003; 
Gonzalez et al., 2008). Many factors have been suggested to explain 
the occurrence of emerging viruses in bats such as diversity (in 
term of species and ecology), ability to fly, immunology, physiol- 
ogy, seasonal migration and roosting behaviour (Dobson, 2005; 
Kunz and Fenton, 2005; Calisher et al., 2006). Another stated that 
population genetic structure correlates with viral richness (Turm- 
elle and Olival, 2009). In the case of Betacoronavirus, the role of: 
(i) the fragmented biotope observed today (the Thai forest biotope 
has drastically reduced over the past decades, leading to the forma- 
tion of species refuge areas promoting species interactions (Harvey 
et al., 2006; Trisurat et al., 2010)) and, (ii) the segregation behav- 
iour and the mobility patterns of bats (metapopulations, migra- 
tions), should be investigated for their impact on the mosaic 
pattern of Betacoronavirus. 


6. Conclusion 


This study suggests that Hipposideridae have been underesti- 
mated in the study of coronaviruses in the wild. Both an Alphacor- 
onavirus and a Betacoronavirus (more precisely the sub-group b 
that includes SARS-CoV) were detected in Hipposideridae in 
South-East Asia (Thailand). In contrast to what is reported for Rhi- 
nolophidae, Hipposideridae are able to host Betacoronavirus close 
to SARS-CoV over long periods of time, suggesting that they could 
be a source of virus for other species under specific circumstances. 
Moreover, the phylogeny of the RdRp, when compared to the 
biogeography, suggests that these two bat families may share an 
ancient relationship with independent lineage of Betacoronavirus. 
Consequently, and by contrast to what have been generally re- 
ported for coronaviruses, we suggest that the host phylogeny glob- 
ally drive the pathogen clustering in the particular relationship 
between Betacoronavirus and Rhinolophoidea. As for other emerg- 
ing viruses, interactions between wildlife refuges (whether consid- 
ered in spatial or organic dimensions) and human populations, 
likely enhanced by fragmentation, changes, cultural behaviour 
and secular traditions, contribute to drive the species barrier cross- 
ing of coronaviruses to humans. 
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